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We study birational transformations of the projective space originating from lattice statistical 
mechanics, specifically from various chiral Potts models. Associating these models to stable patterns 
and signed-patterns, we give general results which allow us to find all chiral g-state spin-edge Potts 
00 ■ models when the number of states q is a prime or the square of a prime, as well as several g-dependent 

C^) ' family of models. We also prove the absence of monocolor stable signed-pattern with more than four 

states. This demonstrates a conjecture about cyclic Hadamard matrices in a particular case. The 
£N| i birational transformations associated to these lattice spin-edge models show complexity reduction. 

In particular we recover a one-parameter family of integrable transformations, for which we give a 
matrix representation when the parameter has a suitable value. 
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INTRODUCTION AND PRESENTATION 



4h . In a previous publication [1] a set of birational transformations acting on projective spaces of various dimension 
Oh- has been introduced. These transformations are birational realizations of Coxeter groups. They arise naturally in 
^ (— i lattice statistical mechanics in relation with the Yang-Baxter equations for solving vertex models and star-triangle 
relation for solving spin model [2, 3]. However it is important to note that these birational symmetries are actually 
symmetries of the lattice of statistical mechanics models beyond the Yang-Baxter integrable situations: they can 
be seen as (generically infinite) discrete and non-linear symmetries of the parameter space of the model and, for 
instance, of the phase diagram of these lattice models [4]. These transformations can in fact be considered per se, as 
, discrete dynamical system. The degree complexity (or entropy) of these transformations has been evaluated [5-7]. An 
' unexpected complexity reduction has been found, and it has been conjectured that the most general Potts model has 
s * ! the same complexity as the most general cyclic Potts model. Considering that the most general cyclic Chiral Potts 
model has only q homogeneous parameters while the most general Potts model has q 2 homogeneous parameters the 
equality between their respective complexities is not obvious. In this paper we go further and study many particular 
cyclic spin-edge Potts Models and their associated birational transformations. As it is described below finding spin- 
' edge cyclic chiral Potts Models for lattice lattice statistical mechanics amounts to finding the so-called stable patterns. 
This problem turns out to be related to many interesting field of mathematics: Bose-Mesner algebra [9], association 
schemes [8], Hadamard matrices (one of our result demonstrates a particular case of a conjecture about Hadamard 
matrices [10, 12]), Gauss identity, etc . 

A y-state Potts model [13] is completely defined by a lattice and a Boltzmann weight matrix W. The spins, which 
have q states, are located at the vertices of a graph, with oriented edges. The Boltzmann weight of a given spin 
configuration is the product of the Boltzmann weight over all edges, hence the name spin-edge model. The Boltzmann 
weight of the edges can be conveniently seen as a matrix : the rows refer to the beginning i of the oriented edge, and 
the column to the end j. The weight of the oriented edge is W(<Ti,<Tj). The so-called inverse-relation [18, 19] 
implies a functional relation between the partition function of a model associated to a matrix and the model associated 
to its inverse for the same lattice. 

By definition a chiral Potts model is a model for which the entries W%j and Wji are different, i.e. the Boltzmann 
weight matrix is not symmetric. Of particular interest are the cyclic chiral models for which the Boltzmann weights 
Wij are functions of i — j mod q. This class contains in particular the integrable chiral Potts models [14, 15]. The 
global symmetries of the cyclic models have been classified in [16] . The most general cyclic Potts model correspond to 
the case where there is no other constraint. It means that the Boltzmann weight matrix W is cyclic. Now we can look 
for other less general models, obtained by imposing further constraints on the entries of W. The simplest constraints 
are equalities between some entries of the matrix W, but we will also consider the case where these constraints are 
"anti-equalities", i.e. we demand that some pairs of entries are opposite. Imposing that two Boltzmann weights are 
opposite could appear unphysical but, as it will be explained below, it is mathematically very natural. 

However these constraints need to be compatible with the inversion relation mentioned above, as well as with 
a Hadamard inverse described below. The aim of this paper is to find such matrices and the associated birational 
transformations. It is organized as follows: we first recall some definitions and what is already known on this problem. 
We then generalize the notions used in this framework and gather together the notations we use. The next two sections 
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are devoted to our analytical results. These results are grouped into two sections, depending whether the result is 
directly of importance from the lattice statistical mechanics point of view, or not. In sectionll we first present our more 
mathematical results. In contrast, the results corresponding to lattice statistical physics are all given in the subsection 
III A while the proofs, together with comments and examples, are given in subsection IIIB. The content of the next 
section III is more "Potts model oriented" and as far as (Boltzmann weight spin-edge) matrices are concerned, focused 
on particular subcases of cyclic matrices, with an attempt to perform an exhaustive classification of the "interesting 
spin-edge Potts models". We give all chiral Potts models when the number of colors q is a square or the square of a 
prime. In some cases we also study the degree-complexity of the birational transformations canonically associated with 
these subcases of cyclic matrices. Here also results are gathered in the subsectionlll A and the proofs in subsection 
IIIB. Then, in section IV, we turn to the cases where we were not able to find analytical results and introduce a 
computer-aided method which enable us to perform some calculation despite the huge combinatorial of this problem. 
We then present these results. 

As far as applications to spin-edge g-state Potts models are concerned our main results correspond to finding the 
stable patterns, or in other words, the "interesting lattice statistical mechanics spin-edge models", when, q, the number 
of states of the g-state Potts model is a prime, or the square of a prime, and providing some first steps for results when 
the number of states is the product of two primes. We also have, as a byproduct, other more specific results, like a 
demonstration of a conjecture about Hadamard matrices (for monocolor stable patterns). Beside the rigorous proofs, 
we give numerous examples, most of them in the appendices. We have tried to be rigorous in the demonstration but 
pedagogical in the examples. The reader more interested in the "lattice statistical mechanics point of view" can skip 
the more mathematical section II, in particular the lemma which are more technical. However we feel that the results 
of this section are worthwhile per se, and can be usefull to go further in the classification of the Potts models. 

A. Recalls 

1. The context 

Starting from the lattice statistical mechanics point of view, we consider an anisotropic Potts model on a square 
lattice with Boltzmann weight matrix Wh for the horizontal edges and W v for the vertical edges. It has been shown 
[18] that if T(Wh, W v ) is the transfer matrix of this model then 

T(W h ,W v )T(W-\j(W v )) = C(W V )I 

where J(W) designates the matrix which entries are the inverse of the entries of W (see below). Transporting this 
equality to the eigenvalues of T permits to find a functional relation for the partition function of the model. This 
functional relation induces a constraining symmetry for the phase diagram. 

We now adopt a more general point of view and we consider the q x q matrices projectively as elements of <CP q 2_ 1 
(since Boltzmann weight matrices are defined up to a multiplicative constant). Using the same notation as in ref [20], 
we define K = I o J where I is the usual matrix inverse I(M) = M^ 1 and J is the Hadamard inverse (inverse of the 
Hadamard product) defined by (J(M))^ = -p— . The transformations I and J are two non commuting involutions, 
which can be represented polynomially in CP g 2_ 1 . In this representation I replaces each entry of M by its cofactor, 
and J replaces each entry by the product of all other entries. It is clear that K and its inverse K^ 1 = J o I are both 
rational transformations. At each step, the q 2 entries of the matrix M are factorized as products of polynomial with 
integer coefficients, and the common factors of all the entries are discarded. 

2. Degree complexity 

A quantity characterizing the complexity is the degree complexity A [5-7] . We simply recall the definition 

A = lim — log d„ 

n — >oo Ji 

where d n is the degree of the n th iterate K n , where K is represented as q homogeneous polynomials of degree d. 
Without the factorizations d n = d n and consequently A = logd. For some transformations one has d n ~ 5 n with 
S < d, this is called a complexity reduction [23, 24]. When the growth of the degree is polynomial one has A = and 
the transformation is integrable[17]. Finally we define the degree generating function as 



oo 

f( x ) = X! dnX 

n=0 



(1) 
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when this serie has a positive radius of convergence. 

3. (I,J)-stable patterns 

We consider a disjoint partition P = {Eq,--- ,E r _i} of the indices with |_l[-=o = {(hj): h3 — 0, • ■ ■ ,n— 1} 
where the symbol |J denotes the disjoint union, we call r the number of colors, and we consider a matrix M such that 

(i, j) e E k and (i',f) G £ fc => M« = M VJ . 

A matrix verifying this set of equalities is said to belong to the pattern P as in ref [1]. We are interested in the 
matrices belonging to the same pattern as their image by K. Therefore we will consider pattern containing at least 
one invertible matrix (for example we exclude the pattern where all entries are equal). Obviously the transformation 
J is compatible with any pattern. Therefore a matrix and its if —image will belong to the same pattern iff this 
matrix and its inverse belong to the same pattern. Such a pattern is said inverse-stable. The number r of subset of 
the partition is called the number of colors. The transformation associated with a stable pattern with r colors acts 
on CP r _i. 

4- Cyclic matrices 

Consider the set of q x q cyclic matrix M c with entries M c (i,j) such that 

M c (i,j) = M c (0,i-j mod g) (2) 

The corresponding model of lattice statistical mechanics is the cyclic chiral Potts model[14, 15]. Chiral refers to the 
fact that M c (i,j) is not necessarily equal to M c (J,i) (the lattice has an orientation) and cyclic refers to the fact 
that one restricts to cyclic Boltzmann weight matrices. The corresponding pattern is inverse-stable and also matrix- 
product stable (see next sub section). Since a cyclic matrix is fully determined by its first row, the transformation K 
can be represented in CP q _i. It is found that complexity reduction does occur for cyclic matrices and the complexity 
is the largest root of x 2 + (2 — (q — 2) 2 )x + 1 [21]. From numerical analysis it has been conjectured that this value 
for the algebraic complexity is the same than for arbitrary matrices without any constraint on the entries. Another 
inverse-stable pattern is provided by cyclic and symmetric matrices. In that case the complexity reduction is even 
bigger and the complexity is the root of largest modulus of x 2 + (2 — (p — l) 2 ) x + 1 where p = [II + 1 where |_ 
denotes the integer part. One aim of this paper is to find some subspaces where further complexity reduction takes 
place. 

B. Generalizations 

1. Product- stability 

A pattern is said product-stable if the product of two matrices belonging to this pattern also belongs to this pattern. 
Using the Cayley-Hamilton theorem, one can express the inverse of a matrix M as a linear combination of its q— 1 first 
powers. Therefore product-stability implies inverse-stability. We are going to present examples where the reciprocal 
proposition is wrong. From now on we call P-stable a pattern which is product-stable, I -stable a pattern which is 
inverse-stable, and /P-stable a pattern inverse-stable but not product-stable. An obvious example are the symmetric 
matrices which are inverse stable (the inverse of a symmetric matrix is symmetric) but are not product-stable (the 
matrix-product of two symmetric matrices is not necessarily symmetric) . 

2. Generalization of the notion of pattern : signed-patterns 
We generalize the notion of pattern and look for a set of r independent q x q matrices Mi such that 

(r-1 \ r-1 

t=0 / i=0 
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Let us introduce the characteristic function x which associates to each set of indices E the matrix x(E) defined by 
{x(E))ij = | g ^ £ ' ^ G P atterns defined in the previous paragraph correspond to 

M k = x(E k ) (3) 
For reasons explained below, we also consider the more general case of matrices with entries 0, 1 or -1 

M k = X (E+)- X (E^) (4) 

We call the partition {Eq,Eq, • • ■ , E^_ 1 , E~_^ a signed-pattern, and r the number of colors. The algebra generated 
by these matrices is J-stable. The notion of "signed-patterns" simply corresponds to the notion of patterns defined 
by equalities between entries up to a sign. 



3. Stability of signed-patterns 

A product-stable set of matrices with entries or 1 and which sum up to the all-one entry matrix is called an 
association scheme [8]. It is an algebra. If the matrices are also symmetric, then it is a Bose-Mesner algebra [9]. 
The problem we address in this paper can be summarized as finding r q x q matrices Mi with entry 0,1,-1 verifying 



1 Vj,fc 



and 



£ y* M > = J2 z > M * ( 5 ) 

\i=l / \j=l / i=l 

for P-stability and 

(r \ _1 r 

^XiMA =J2ziMi (6) 
i=l ) i=l 

when the inverse of Xw=i x i^'h exists for /-stability. From the definition, the matrices induced by a F-stable pattern 
form an algebra. But the matrices induced by /-stable patterns do not always (in contrast with the problem studied 
in [16]). 

Note that if a set of matrices {M,} defines an algebra, so does the set {P~ 1 M i P (T } where a is a permutation of 
{0, • • • , q — 1} and P the associated permutation matrix (Pa)^ = ^i,a(j)- However if M is a cyclic matrix, P~ X MP„ 
is not necessarily a cyclic matrix. 



C. Notations and definitions 



1. From now on we will restrict ourself to cyclic matrices. As far as notations are concerned, we identify a cyclic 
matrix and its first row seen as a vector in CP q _i. Let v G C q be a vector, we use the notation Cy(v) to denote 
the q x q matrix 

(Cy(v)) y . = Vi-j 

and Diag(v) to denote the q x q matrix 

(Diag(v)) y = ViSij 

6 is a Kronecker symbol. 
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2. As already mentioned in the introduction, the discrete Fourier transform plays a crucial role for stability of 
cyclic matrices. We therefore define the matrix U = (w 1 - 7 ) where u = cxp ^-i and use notation 

x = Ux 

to denote the Fourier transform of the vector xeC. We note the relation 

U* x Cy(x) xU = q Diag(x) (7) 

which will be useful. For the reader familiar to lattice statistical mechanics U is the generalization of the 
Kramers- Wannier duality, however it is not a transformation of order two, but of order four. 

3. When patterns are explicitly given we use a straightforward representation: we put in a bracket the entries of 
the first row. An example is given with the comments in appendix A. 

4. The subspace spanned by the set of vectors {v(i)} i=1 ... r with complex coefficients is denoted: 

Cv(z) 

l<i<r 

5. We use arithmetic modulo q, Z* is the set of elements of Z 9 which are invertible and if d is a divisor of q (noted 
d | q) one introduces 

Z(q, d) = {k G Z q | gcd(fc, q) = d} 

6. We will also use the convolution product (noted *) and the Hadamard product (noted .) between two vectors 
of C 9 defined respectively by 

9-1 

( U * V )* = ^ U 3 V i-j 
3=0 

(u • V) ; = UiVi 

It is straightforward to see that Cy(u)Cy(v) = Cy(u*v), With these notations, Eq. 7 reads: 

u7v = qu ■ v 

We note u*™ the convolution product of u with itself n times. By convention u*° = X ({0}). Keeping in mind 
that diagonalization of the cyclic matrices requires the discrete Fourier transform Eq. 7, the previous relation 
amounts to writing that the eigenvalues of the product of two cyclic matrices is the product of the eigenvalues. 
With these notations a partition £ — {E\, • • • , E r } is P-stable if Vo*, h, 3a such that 

OiXiEi) * J2 b iXi E i) = CiX&i) 

i i i 

the partition £ = {E± 5 • • • , £^ r }is /-stable if Vai, 3ci such that 

a lX {E t ) * J2 ^Ei) = (1, 0, • • • , 0) 

i i 

corresponding to the matrix inversion of a cyclic matrix. 

7. A set of disjoint subsets £ = {E = {0} , E x , ■ ■ ■ , Ek} of {0, • • • , q — 1} is convenient if V (m, ■ ■ ■ , rife) G N fc V/ G 
[O.fc] 



V*, j G E l (x(Bi)*" 1 * • • • * x{E k )* n % = (x(^i)*" 1 * • • • * x(E k ] 



*n k \ 

h 



By a slight abuse of notation a set E such that {{0} , E} is convenient is also called convenient. In that case 
one has 

Vn>0 Vi,jeE (Cy(x(E)r) 0ti = (Cy(x(E)r) 0J (8) 

Actually, since any power of a q x q M can be expressed as linear combination of the first q — 1 powers with the 
help of the Cayley-Hamilton theorem, one needs to verify Eq.8 only for < n < q — 1. 

Intuitively, a convenient set of disjoint subsets can be seen as a possible "beginning" of a stable pattern. Indeed 
it verifies some necessary conditions such that it can be it extended to a stable pattern. In particular each set of 
a stable partition is convenient. This will be used in section IV. Note that if £ is a partition then it is P-stable. 
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8. An admissible set is a subsei E of {0, • • ■ ,17 — 1} 

E = |J d i d H d (9) 

where Z? is a subset of the divisors of q, is a subgroup of Z\ and gcd(did, q) = dVd e D. We will show below 
that the union in Eq. 9 is indeed a disjoint union and that the admissible sets are in fact the only possible sets 
in a stable pattern. 

9. For I a finite set, x = (xi)i e i we note V x j = {A±, ■ ■ • , A r } the partition of I such that Xi = Xj iff i and j are 
in a same A^. Actually the A^s are the preimages of the application i — > Xi. When the set I = {1, • ■ • , q} we 
do not specify it and note simply V x . For example for q = 4 Vr XllXl ,xi,xi) = {{1, 3, 4} , {2}} 

10. Finally if E is a group, F < E means that F is a subgroup of E. 

II. ANALYTICAL RESULTS PERTAINING TO MATHEMATICS 

In this paragraph we express the inverse-stability and the product-stability for the pattern and the signed-pattern 
in term of matrix subalgebra and list mathematical results which are used in the next section. We also present results 
interesting per se, and likely to be usefull to go further in the classification of the lattice models. 

A. List of the results 

1. Pattern stability as matrix subalgebra 

• The pattern £ = {Ei} is product-stable iff there exists a partition T = {Fj} such that 

C^E~) = C X (Fj) (10) 

l<i<r 

• The pattern £ = {Ei} is inverse-stable iff there exists a partition T = {Fj~ , F~] such that 

c( X (F+)-x(F7)) (11) 

l<z<r 

• The signed-pattern £ = {E^~ , is product-stable iff there exists a partition T = {Fi} such that 

C fe?) - x(if)) = Cx(^) (12) 

l<i<r ^ ' l<j<r 

• The signed-pattern £ = {Ef,E~^\ is inverse-stable iff there exists a partition T = {F^ , F~} such that 

C fe+j - x5f)) = C ( X (F+) - X (F7)) (13) 

l<i<r ^ ' l<j<r 

2. Stability by multiplication 

If £ and T are two signed-patterns verifying Eq 

a£ ■ 

By a£ = £ we mean Vi, 3fe such that either aEf - 



. 13 then for any a prime with q 
= £ and aT = T 

= £+and aEl = E~ , or a£f = EV and oJST" = E+ . 
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3. Admissible subsets 
All the sets E i or Ei in relations Eq.10-Eq.13 are admissible. 

4- Convenient sets : necessary conditions of stability 
This result is mainly useful for the demonstration of the result of the paragraph III A 2. 

For a partition {E\, ■ ■ ■ , E r , A} we define a partition {Fi, • • • , F s } such that ©i<j< r Cx(-Ej) C 1<3<8 Cx(-Pj), 
with Fj maximal. We define J the subset {1, • • • , s} by 

j g J ^ fc ^ o Vz g {i, • • • , r} , fc e F s 

Then the set {E 1 , ■ ■ ■ , E r } is convenient i/f 

0C^ C I CxCEj) ] f0Cx({a})) 

j£.I \l<i<r I \a£A / 



5. Subgroups induce product-stable patterns 

If H is a subgroup of Z* (the set of the invertible elements of Z g ) then H induces a product-stable pattern given 
by the classes modulo H (i.e. the{iH} i e Z q ) 

6. Monocolor inverse-stable signed-pattern 

Except a q = 4 (and q = 1) example, there is no inverse-stable signed-pattern. This result proves a conjecture 
concerning cyclic Hadamard matrices in the particular case of symmetric matrices. 

B. Proof and illustration of the above result. 

Below we prove and comment the results mentioned above. We also give examples and illustrations. First we need 
two assertions. 

Assertion 1: Let V C C q be a vector subspace of dimension r then V is product- stable iff there exists disjoints 
subsets F u ■ ■ ■ , F r of {1, • • • , q) with V = ©!<,<,, C X (Fi). 

Assertion 2: Let VcC be a vector subspace of dimension r and V* the (supposed nonempty) subset of vectors of 

V with all non-zero components, then (V*) C V iff there exist a partition of {1, • • • , q} F^ , F^ , • • • , F+, F~ 
with 

V= C(x(F+)- X (Fr)) 

l<i<r 

In other words Assertion 1 states that V can be generated by vectors with entries or 1, and Assertion 2 states 
that V can be generated by vectors with entries 0, 1 or -1, and such that the absolute value of all the entries of these 
vectors sum up to the all one entry vector. These assertions are proved by recurrence on the space dimension q. 

Proof of Assertion 1 For q = 1 the assertion is clear. Suppose Assertion 1 is true for dimension q — 1 and let 

V e C q be a product-stable subspace. Let V fl (C^ 1 x {0}) = W x {0} therefore W = ®i<;< s Cx(i ; ]) where the F t 
are disjoint subsets of {1, • • • ,q— 1}. Let v £ V, we will show that i, j £ Fi implies m = Vj. If v q = it is clear. If 
not, one can always admit v q = 1, then v • v — v G W x {0} so that — Dj = — Vj, which implies i>j = vj since 
Vi + vj — 1 is impossible (one can add xffi) x {0} to v). So if v G V is a vector such that v q = 1, we can admit v i = 
for i G Fi VI. Since V = (W x {0}) © Cv and v v = v then all v = x(F) where F C {!,-■■ ,q} disjoint of all Ft . 
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Proof of Assertion 2 Let V satisfying the conditions of Assertion 2 and v G V* such that v q = 1, we define 
W = {weV|u., = 0} and i" = {« | to* = Vw G W}. If J = {1, • ■ ■ , g} then W 7 = Cv therefore v" 1 = v which 
proves the assertion. On the contrary if I ^ {1, • " > there exists w G W such that u>; ^ for all i £ I, this is 
possible while W is not a finite union of proper subspaces. Let u = lim e _> (v + e _1 w) , m = v7 if i € J and 
Ui = if i ^ 7. On another hand v — v _1 € W so = w~" for i G 7, finally Vj = ±1 for i e I. One has V = Cu © W, 
and we proceed again with W, keeping only the coordinates not in I. 



1. Proof of II A 1: pattern stability as matrix subalgebra 



Eq. 10 and Eq. 12 are direct consequences of the assertion 1, taking 

l<i<r 

Eq. 11 and Eq. 13 are direct consequence the assertion 2. If Eq. 10 holds, then {0} G £ and {0} G P. Indeed if k ^ 
and are in the same Fj and 1 G -Ej then ^x(Ei)J = (x{Eifj is impossible. 

We give in Appendix A the exhaustive list of theP-stable, /-stable and IP-stable pattern for q = 8. The pattern 9 
is an illustration of 10, the pattern 5 is an illustration of 12, the pattern 18 is an illustration of 11, and finally pattern 
8 is an illustration of 13. 

Note that inverse stability is the justification of introducing signed-pattern, which could not be justified in the 
strict framework of lattice statistical mechanics since asking that two Boltzmann weights are opposite is unphysical 
(however such opposite entries in the Boltzman weight matrix occurred in the solution of the 3d generalization of the 
Yang-Baxter equation (tetrahedron equations) by Baxter and Zamolochikov) . 

In Appendix A we give a detailed non trivial example of application of Eq. 13. 



2. Proof of II A 2: stability by multiplication 

Let £ and P be two signed-patterns verifying Eq. 13 and A = ©i<i< r C ^x(E^) — x(-E^ )J = 
0i<j<r C (x(Fj) - X{FJ)) therefore for any i, ^E+) - X (£7) g A. This implies that, for any j, (x(ifj) fc = 

e ^X(E i for k, I £ Fj~ U FT with e = 1 if k and I are both in F^ or both in FT , and e = — 1 else. Let us introduce 
the polynomial 



G 



p{x) = J2 xek - £ xek ) - e £ xel - xe 

One has P{u>) — and consequently, using a Galois symmetry argument, P{uj a ) = for a prime with q. Therefore 
x{aEf) — x(aE7) G A and © 1<i<r C ^x(E^) — x(E~)*J C A. The equality follows by a dimension argument. 
Notice that, applying this result to a = —1, one finds that if E G £ then either £7 = — E or — £7 is another set of £ . 



3. Proof of II A 3: admissible subsets 



Let us first recall that an admissible set E of {0, • • • ,q— 1} is a disjoint union 



E = |J d i d H d 

where D is a subset of the divisors of q, Hd is a subgroup of Z% and gcd(di £ j, q) = dVd G D. 

We first note that the intersection of two admissible sets is an admissible set. This comes from the fact that if d 
is a divisor of q, H and H' are two subgroups of Z\ , and id and i'd are two elements of Z(q, d) then either there 
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exists i"d e idH n i'dH' implying that idH n i'dH' = i"d (H n H'), or idff n i'dH' = 0. We need also the following 
technical lemma: 

Lemma 1 Let P(X) G Z [X], »/ id G Z(q, d) then {jd G Z(g, d) | P(w id ) = P(ui jd )} = idH with H < Z% (with a little 
abuse of notation). 

Proof: uo ld is a |th primitive root of unity that we denote £. Let k be the inverse of i modulo \. The condition 

P(uj ld ) = P(u jd ) becomes P(() = P(C kj ). But the set H = jm G Z% \ P(Q = P{C m )} is a subgroup of Z|. Indeed 

P(C) = P(C m ) does not depend of the particular choice of the primitive root of unity since it amounts to saying that 
P(X) — P(X m ) is a multiple of the minimal polynomial of C and therefore P(C) = P(( m ) and P(() = P(( l ) implies 

P(0 =P(C ml ) andP(C) =P(C m "). 

Let us consider Eq. 13. Let k G P+ and define P l (X) = E eeS + x& - E e e£- X<i • B y Ec l- 13 one has 

P+= f| {iGZ fl |P i ( W ')=i'i(w fc )} 

l<i<r 

if d = gcd(o,fc) then Ff R Z(?, d) = fli<i<r {* e Z(g, d) I Pi(w') = Pi(^ fe )} which is an admissible set by the results 
above, and finally Fj~ is also admissible. 

4- Proof of II A 4: convenient sets, necessary conditions of stability 

We define the intersection of two partitions £ = {E\, ■ ■ ■ ,E r } and T = {Pi, • • • , F s } of a finite set I by 

£ nF= {Ei n Fj |1 < i < r 1 < j < s} 

Lemma 2 Let y £ C" and a(l), ■ • ■ , a(t) G C" /or A G Pa(i)H ■ • • PI P a (t)- 27ie iwo following affirmations 14 and 15 
are equivalent (with the convention — 1) 

n 

V(fei,--- ,fc t )€N 4 ^a(l)* I ..-a(t) 1 fc *y i = (14) 
i=i 

VA g P a( i) n • • • n p a( t) »* = ( 15 ) 

Proof: The proof goes by induction over f. For 4=1, let us define the sets A, by P a (i),{i, -- ,n} = {^li • " ■ j A r } and 
let 6j be the value of a(l) 4 for i in A,-. Using 14 one has 

(■- 1 

p 

Since detP = rii<i<j<r(^i ~ fy) ^ (Vandermonde determinant), EigA fc Vi = f° r any 1 < k < r. This proves the 
property for t = 1. 

Let us take £ > 1 and assume the lemma for t — 1, therefore Eq. 14 is equivalent to 

X)a(*)h/i = VAGP a(1) n---np a(t _ 1) Vfc t 6N 

and using the lemma for t = 1 this is equivalent to 

5^ W = VC e ^a(t),A With A G Pa(l) (1 • • • Pa(t-l) 



\ Ez G A r w / V 



(16) 
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^z/i = o CeP a( i)n-nP a(t) 

ieC 

which completes the proof of the lemma. 

We now use this lemma to prove the result II A 4. Let {Eq,E\,--- ,E r ,A} be a partition such that £ = 

{Eq,Ei,--- ,E r } is convenient. For (fcor" ,k r ) e W +1 one introduces u = x(Eo)* k ° • • • x(E r )* kr ■ By definition 

Ui = Uj for i and j in the same E%. If v = u = q k "^ ^ kr x(Eo) ■ ■ ■ x(£V) this implies Yl m "mW" ,m = J2 m v m Lj~ jm 
leading to 

£ (u- im - w-' m ) (xlEo)) k ° ■ ■ ■ (xm) k " = 

— \ / m \ / ni 

m 

We now use the lemma 2 with y m = uj~ lm — uj^^ m and a(n) = x(E n ) to get the result. 

5. Proof of II A 5: subgroups induce product-stable patterns 

To show that subgroups induce product-stable pattern, we first note that the class modulo H is indeed a pattern. 
Let I be a set of representative of the classes modulo H. Let 

x = £a,x(iff)e0Cx(iff) 

iei iei 

then 

iei ' ' heH 

where \A\ denotes the cardinality of the set A. Consequently if j' = tj with t £ H then (x)^ = (x) jV which implies 

0C^)C0C X (^) 

iei iei 

The inverse inclusion is shown using inverse Fourier transform. 

6. Proof of II A 6: monocolor inverse-stable signed-pattern 

Let us take q > 1 (the case q = 1 is obvious). It is readily verified that the monocolor signed-pattern [a, —a, —a, —a] 
is /-stable. We now prove the following lemma which we will need. 

Lemma 3 i/o,i*eN* o-^ab = ^K^V we means the elements of 1 a b which are prime with b). 

Proof: Let us introduce the set C = {k G 2, a b \ gcd(fc, b) = 1}, we will prove aZ* b = aC. It is clear that al* ab C aC. 
We need to show that if k is prime with b, then one of the k, k + b, ■ ■ ■ , k + (a — 1)6 is prime with a. Let us write 
a = cd where the prime factors of b appear only in c. Since gcd(6, d) = 1 then k, k + b,---, k + (d — 1)6 are distinct 
modulo d, therefore one of them is equal to 1 modulo d, which proves the lemma. 

Below we show that there is no other /-stable signed-pattern than [a, —a, —a, —a]. Let E + , E~ be an /-stable 
monocolor signed-pattern and M = Cy (x{E + ) — x(E~)), the inverse-stability can be expressed as M 2 = tl q where t 
is some even non zero integer and I q is the identity matrix. Applying twice M to the all-one entry vector, one gets 

M 2 = s 2 I q (17) 

with s = \E + \ — \E~\ where we consider, without loss of generality, that \E + \ > \E~\. 
We now prove that s 2 = q. Indeed using Eq. 13 

U*MU = q diag (£(E+) - xjE- 
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so there exists a constant k and a partition {F + , F } of {0, q — 1} such that x{E + )—x{E ) — k {x{F + ) — x{F ))■ 
s = (x(E+) - X(i r )) = k (x(F+) - x(F-)) so a = k 

xiE+) - xl^n = a ( X (F + ) ~ X(F~)) (18) 
Define N = Cy (x(E + ) — x(F~)) applying again the Fourier transform to Eq. 18 one gets 

which combined with Eq. 17 yields 

M 2 = ql q 

as stated before. 

Applying the equation above to a diagonal term proves that M is symmetric and therefore Eq. 17 can be written 

as 

MM = ql q (19) 

In other words, M is a so-called [12] symmetric Hadamard matrix, see ref [11]. 

Using Eq. 17 and Eq. 19 one gets q = 4u 2 . The example given in the beginning of this paragraph corresponds to 
u = 1. 

Since X (E+) - x{E~) = 2 X (E+) - x ({0, 1)}) = 2 X (E+) - q X ({0)}) and using Eq. 18 



K X(E+)) = ±u for i ^ (20) 
Since E + is an admissible set (see paragraph II A 3) 

e + = y di d H d 



deD 



where D is a subset of the divisors of q, gcd(d i d ,q) = d and H d < Z% . Using the result of subsection II A 2 and since 

\E + \ 7^ \E~\ for any a £ Z* one has aE + = E + yielding 1>* q E + = E + , in particular for d G D, i d dZ* C E + . Using 
Lemma 3 one has i d dZ* = i d d7h\ = dZ* q , so that 



E+ = |J dZ| (21) 



deD 

(in the example shown in the beginning of this section one has E + = {1,3} U {2}) . 

At this point we need to use results of number theory. The so-called Moebius function \i [29] is defined by 

1 if 71 = 1 

M n ) = { ( — 1)' if n = pi ■■■ pi with pi, ■ ■ ■ ,pi distinct primes 
else 

and it has the property that X^fegz* C k = M( rl ) where C = exp — i. In our case one gets X^fcgz* Ljjkd = /^( § ) • We now 

n q/d 

use Eq. 20 and Eq. 21 and we get 



±u=(x(E+)) i = J2»( 1 d ) 



deD 

„2a, 2a, 



Noting q = 2 2a °p 1 ai ■ ■ -p^ 1 where 2, p 1: ■ ■ ■ , p l are distinct prime numbers, one has 

\{t\q such that /x(t) = 1}| = 2 l = \{t\q such that fj,(t) = 
therefore u = 2 a °~ 1 p a l 1 ■ ■ - pf l <2\ consequently I = and Oq = 1, which proves that q = 4 is the only possible value. 
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III. ANALYTICAL RESULTS PERTAINING TO STATISTICAL MECHANICS 



In this paragraph we express the inverse-stability and the product-stability for the pattern and the signed-pattern 
in term of matrix subalgebra and list our analytical results. 



A. List of the results 



1. q is a prime 



When q is a prime number, there is no other product-stable pattern than the patterns induced by the subgroups of 
Zq. Furthermore there is no inverse-stable pattern which is not product-stable. Consequently there are 1 + r(q — 1) 
stable patterns, where r(n) is the number of divisors of n. 



2. q is the square of a prime or the product of two primes 



If q is the square of a prime p, then there are 1 + r(p — 1) + r 2 (p — 1) product-stable patterns. In addition to the 
class modulo Z q described in II A 5, there exists other product-stable patterns. If P denotes the natural projection 
from Jj q on Z p the set of product-stable patterns is: 

{{0} , Z, \ {0}} (J {aL} (J {{0} , aP- 1 (H), bpK} 

where 



a£Z, L < 1j q 

a,b£Z* H, K < Z* 



If q = pip 2 is the product of two different prime numbers p\ and p 2 then the three-color pattern defined by Eq = {0} 
and Ei = {pi,2pi, ■ ■ ■ ,(p 2 — is product-stable. 



3. An integrable one-parameter family of integrable patterns 



if q 2 is a prime the three-color pattern formed by 

E a = {0} 

Ei = {i 2 ,ieZ* q } (22) 

E-2 = Z q — E() — El 

is product-stable. Furthermore if q = 4k + 1, then the associated transformation K is integrable. 



4- Six families of stable patterns 



We give below six families of patterns which are stable for any even q. Two patterns on the same row of the table 
below are related by discrete Fourier transform. 

Prod. -stable signed-pattern Inv. -stable simple-pattern 

Pi = [x ,Xi, ■ ■ ■ ,-Xi, ■■ ■ ,Xi] Ql = [x , ■ ■ ■ ,Xq,Xi,Xq, ■■■ ,x ] 
P 2 = [x ,Xi,-Xi, - ■ ■ ,X 2 ,~Xi,--- ,Xi] Q-2 = [x ,Xi,X , ■ ■ ■ ,X ,X2,X ,Xi, ■ ■ ■ ,Xi] 



Pa 



xo,xi,x 2 , ■ ■ ■ ,Xg,—Xl, 



X ,Xi,Xo,X2,X ,X 3 , 



,X ,X9 



Patterns Pi and Qi are two-color pattern, P 2 and Q 2 are three-color pattern , and P3 and Q3 are | + 1-color patterns. 
For pattern Pi (resp. Qi) the entry— x\ (resp. Xi) is in position | + 1 (position starting at zero). For pattern P 2 and 
Q 2 the entry x 2 is also in position § + 1. Note that for P 2 the elements before and after x 2 are a^when | is even and 
—Xi when § is odd. 



13 



B. Proof and illustration of the above results. 

1. Proof of III A 1: q is a prime 
The key point of this demonstration is the well known fact that (if a , ■ ■ ■ , a 9 _x G Q) 

9-1 

^ aiUJ 1 = a = ■ ■ ■ = a g _i (23) 

i=0 

Let £ = {E^, Ff , ■ ■ ■ , E+ , E~ } and T = { F x + , FT , • • • , F+ , F ? 7 } be two partitions verifying Eq. 13 with Ef , ■ ■ ■ , E+ , 
F^, • • • , F+non empty. We admit r > 2 since the case r = 1 occurs only for g = 1 or q = 4 as shown in section II A 6. 
Note that the E's and the F's play the same role and can be interchanged. 
Let 1 < l,m < r, Eq. 13 can be rewritten as 



X (E+) - X (E l ) = X (E+) - X (E] ) ) if i, j G F+ or i, j G F„ 



X (F+) - x(F ; ) = - x(F+) - X(^D if F+ and j G F, 



(24) 



If G F+we will show Fj + = {0} and Fj~ = 0. Suppose, a& absurdo, that there exists i 7^ with i G F+ U F[~ . 
Using Eq. 24 one gets 



X (F+) - X {EJ) ] = ± ( X (F+) - x(FD 



for J ^ 1 and G F+. So that 



l^ + |-l^rl = ± E E wi 

which is impossible by Eq. 23, since ^ iE^~ and «F ; + n iF ; ~ = 0. The same applies interchanging the F's and F's. 
We now show that Ey = ■ ■ ■ = E~ = Ff = • • • = F~ = 0. Recalling that E^ = {0} and Ff = 0, let suppose that 

Ff is non empty and i G F ; + , j G F ; ~, using Eq. 24 one has (x(Ei)j = — (x(Ei)j and consequently 1 = — 1, a 

contradiction. 

So, when q prime, Eq. 13 reduces to Eq. 10 which we write 

(xC$j). = (xSF)) e F, VfcJ (25) 

if 1 £ G £ and 1 G A' G F so Y^heH^' 1 = Y.heH^' 1 f° r e ^ • Using Eq. 23, one deduces iF = jH. Taking 
j = 1 one gets i E H and therefore K C F. Interchanging F and A one finds that K = H . From ij' -1 G H we 
deduce that the subset of the partition which contains 1 is a subgroup. 

Again using Eq. 25 with I > 2 one has iEJ = jEJ for i,j G F and therefore F ; + = FF ; + and EJ is an union of 
classes modulo H. On another hand again using Eq. 25 



(x(F)) . = (x(F)) Vi,jGF+ 



therefore iF = j'F and F ; + is one class modulo H. This completes the proof (noting EJ = {0} = OF). 



2. Proof of III A 2: q is the square of a prime 

We consider the case q = p 2 with p a prime and note u = exp ^-1 and C = exp =^-i. We recall that the cyclotomic 
polynomial of respective order p and q are 

$ p (A) = 1 + --- + AP- 1 
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therefore if P(X) = J2ll Q a kX k and Q(X) = £Lo b k* k are two polynomials in Q[X] then 

P(C) = <s=> a = ■ ■ ■ = a p _i (26) 
Q(w) = (p | (/ - m ) 6, = V) (27) 

A minimal polynomial of lu over Q[£] is X p — £. Finally we note P the natural projection Z g — > Z p . 

Lemma 4 If H < Z* then either H = P _1 (A) wii/i A' < Z* ('in £/ia£ case P is sairf to be of type 1) or 
h G P, Zi' e P, Zi 7^ hi implies p\{h! — h) (in that case H is said to be of type 2). 

Proof: Suppose H < Z* is not of type 2, then there exists h, h', h' ^ h such that p | (h — hi), so h!h~ x = 1 + kp with 
1 < k < p - 1. Therefore (1 + fcp) z = 1 + kip for all Z and 1 + pZ p C H and finally P(l + pZ p ) = H: P is of type 1. 

The next lemma is devoted to the determination of P^j^. when E = iH (p\i and H < Z*) or A = jpA (p j j and 
A < Z£). 

Lemma 5 One distinguishes the three following cases 

E = jpK, p\j, A < Z* p => V%ft = { P Z P , {iP- x {K)\ p\i}} (28) 

E = iP-\K), K<Z* p , p\it=^V^ ={{$}, l? q ,{pjK\p\j}} (29) 

E = iH, p\i, H < Z*, H of type 2 =^> 7^ = {jP| j G Z,} (30) 

Proo/ of Eq. 28 

If p j i and p \ I then 

'x(P)) = (xC®)) ^ XI - wiiPfe ) = ^ H (C* 7 * - C yfc ) = <^> iA = IK 

1 k£K keK 



on another hand ^x(P)) = \K\ for any t. The equivalence noted refers to 26 with P < Z* p \ i 



Proof of Eq. 29 

If K = Z*, then E = Z* and 29 is verified. We now suppose A ^ Z*. If p { Z then 

(x(P)) ; = E k eK^n= ^ l{mp+k) = °, on another hand if p\t = p£ fceK C* fe ^ using 26 and A ^ ZJ, 

the value X) fee if C* fe depends only on A using again 26. Finally (x(E)) = P \K\- 
Proof of Eq. 30 

If p\l andpj j the two quantities (x(E)J = Y^heH^ 11 an d (^C^O^ = Tl,h^H ujllh are equal ijff ijP = UH. 
Indeed let us suppose ^x(P)^j = f%(P)J then Vfc G [l,p — 1] , |ijP H (fc + pZ p )| =0, or 1 since P is of type 2, and 
if t G (ijP n (fc + pZ p )) \ (UH n (k + pZ p )) then, using 27, i + p, • • • ,t+(p-l)p also belong to 
(ijH n (fc + pZ p )) \ (UH n (fc + pZ p )), a contradiction. On another hand (x(E)j can be seen as a polynomial in u 

of degree strictly smaller than p over Q [£] therefore ^x(P)^) . ^ Q [C] ( see above point 4 of the recalls). Finally if p j Z 

and p { j the two quantities (^(P)) = £ fc ff C yh G Q [C] and (^E)) - ^ fteff C'" 1 G Q [Q are equal iff 
ijP(H) = UP(H) using 26. 



Lemma 6 If E is convenient, admissible, E n Z* ^ and P n pZ* 7^ i/ien E = Z q \ {0}. 
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Proof: Let us take a set A verifying the conditions of the lemma, since E is admissible then E = iH U jpK with 
H < Z* ,K < Z*and p \ i and p \ j . If there is A G P^ with A = pM C pZ* using (*) = (x(A)) , 

therefore Yl,m£M ujtmp = EmeM £<;J ' mp i so SmeJW C* m = [Afj which is impossible. Consequently there exists t 
relatively prime to p such that t,pe A G ^xCE) anC ^ S0 = (^(-^)) which reads 



c jkt = E + E 1 ( 31 ) 



E + E 

TiGif fceiC hei? keK 



So J2h£H 10 ^ *Q [C] anc ^ ^ cannot be of type 2 since the degree of u> over Q [Q is p. 7? is then of type 1, 
H = P- l {L) where L < Z*. Using 31, £j. eK ( jkt = pJ2ieL C 1 + \K\ • We use 26 and we note 
A s = \{k G K\p | jkt — s}\ and B s = \{l G L|p | jZ — s}\ for 1 < s < p, one deduces \K\ = —A s + pB s which is 
possible only if A s = B s = 1 Vs, therefore A' = L = Z* and finally A = Z 9 \ {0}. 

Remark 7/ {{0} , A, F} is convenient then MA G ^^^)^^^f) one ^ as (x(-^-)J = ifi,j€E or i, j G A 

and so ifP^p^ = {-^li ' ' ' 7 A r } ifeen A C P s and A C A t for some s and t. 

Proof of the result III A 2 We consider below five possible cases of couples (E,E') such that {{0}A, A'} is con- 
venient. In each case we are going to apply the remark above with a suitable choice of the set A, as well as Lemma 
5. 



1. E = iH, p \ i, H < Z* , E' = jK, p \ j, K < Z*, H and K are both of type 2. Using Eq. 30, V 



X(E) 

{kH\ k G Z q ) and Vgj* = {kH\ k G Z q } so P^ n P^ = {fc(A n A)| k G ZJ. Taking A = i? n A', 
^xCA) = ^ ^1 ' ^ one ^ e ^ s ^ C K^^K) and A' C m(HC\K) for some Z and to which is possible 
only if H = K. 

2. E = iH, p\i, H < Z* , H is of type 2 and E' = jP~ l (L), L < Z%, p\j, P^ = {{0} ,Z* q , {pkL\ p \ k}} . 

Taking A = H G P-r^nP-T^ , jP~ 1 (L) C mi/ for some to which is in contradiction with the fact that if is 

& x(£) x(-F) 

of type 2. 

3. A = iH, p\i, H < Z\ , H is of type 2 and A' = jpA, K < Z*, P^ = {pZ p , {tP^{K)\ p\t}}. Taking A = 
Hr\P~ l {K) which a subgroup of type 2, Vg~^ = {s(H D P' 1 (K))\ s G Z 9 } and therefore H C sfffnr'fA')) 
and pA C n for some s and t. We deduce that pH = pK. 

4. E = ipK, p\i, K < Z* and E' = jpL, p\j, L < Z*, V^ffPg^ = {pZ p , {mP-^K n L)\ p f m}}. Taking 
now A = P" 1 ^ n L), = {{0} , Z*, {pt(K C\L)\p\t}} yielding K = L. 

5. E = iP- x {K), K < Z*, p \ i and £" = jP-\L), L < Z*, p \ j. Taking A = p(K n L) G F— } nP— } = 
{{0},Z;,{p^nl)| P^ = {pZp^mP-^KnL)] p\m}} one gets A' = L. 

We now take £ = {E x , ■ ■ ■ , E r ] and T = \v ' x(a p^ H)) = {{0} , Z*, {prniJ | p f to}} A i; ■ • ■ , A r | verifying Eq. 10. If 

one of the AVs is of the type of Lemma 6, then £ = {{0} , Z q \ {0}}. From now on we consider the case where no Ek 
is of this type. So all the AVs are either iH, or jpK or {0}, with H < Z* p \ i or K < Z* p \ j. It is clear that if 
E, E' (distinct) are in £, then {{0} E, E'} is convenient. Suppose there is a set Ek = iH with H < Z* p \ i with H 
of type 2, then using the three first points above one gets that the E\ are the class modulo H . Using now the last two 
points, we see that VA G £ \ {0} either A = iP _1 (A), p \ i or A = jpL, p \ j. Thus we have shown that the only 
possible P-stable patterns are those occurring in the paragraph III A 2. 

We now verify that these cases are indeed P-stable. Obviously{{0} , Z q \ {0}}and {aL} are P-stable. Let us 
now take £ = {{0} , aP^fAf), bpK} with H,K < Z* and a,b G Z*, V 1 o p- = {{0} , Z*, { P toPT| p\m}} and 

^xibpK) = {p^p 1 {^A _1 A| Note that these two sets are independent of a and b. Showing that £ is product 

stable amounts to showing that it is convenient. The only possible A in T \(p^i(H))^ : 'jdpK) are P m ^ an( ^ ^P _1 (A') 

with p \ to and p \ t. When A = pmH then V-rr, = {pZ p , {mP _1 (iI)| p f to} } while when A = tp- 1 (K), V-rr-. = 

{{0} ,Z*,{pmHK\ p\m}} . In both cases Eq. 10 is verified which completes the proof of III A 2. 
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We are now in position to compute the number of stable patterns. We recall that a cyclic group of order n has r(n) 
subgroups so that Z* has r(p 2 — p) = 2r(p — 1) subgroups, but there are as many subgroups of type 1 as subgroups 

in Z* i.e. r(p — 1), so there are r(p — 1) subgroups of type 2. Finally | {{0} , aP- 1 (H), bpK) \ = (r(p - l)) 2 , yielding 
the stated result. 

In appendix we give as an illustration the example of all stable patterns for q = 25 = 5 2 . 



3. Proof of III A 3: an integrable one-parameter family of integrable patterns 

• Clearly for the pattern defined by Eq 22. 

card(-Bi) = card(i?2) = g ^ 

• If k € i?i then there exists a such that k = a 2 mod q therefore 

kEi = a 2 {i 2 mod q) = E% 
so J2ieE! w 1S independent of k \/k £ E\. Let us define A and A' 



4= A ' = £ 



w k2i 



where fei G E\ and fc 2 £ £2- One has A + A' + 1 = 0. 
Let us recall the Gauss's result [25] 



H ex P 



J'=0 



(l + i)y^ q = mod 4 

2z7rj :i; J ^/g 9 = 1 mod 4 

g = 2 mod 4 

i\/<7 5 = 3 mod 4 



from which one deduces 



with 



A = 



1 if q — 1 mod 4 
z if g = 3 mod 4 

• Finally, the Fourier transform is represented by the 3x3 matrix: 

/I ¥ ¥ 
F, = 1 - e "f- X ] (32) 

When g=l mod 4 then e 9 = 1, Eq. 32 corresponds to Eql7 and Eq 18 of [26]. We deduce that the integrable 
mapping discovered in this reference corresponds to a spin edge model of lattice statistical mechanics when q is a 
prime number with q = 1 mod 4, the pattern is explicitly defined by Eq. 22. The homogeneous expression of the 
transformation K : (x, y, z) — > (X, Y, Z) can then easily be found. If one introduces then the inhomogeneous variables 
u = j and v = - a if— invariant having a particularly simple form is 

(u — v) 2 ( n uv — u — v + V 



q + 2- 

(2uu — u — f )(u + i) — 2) \ ti + u 

When g = — 1 mod 4 then t q = 1 and the corresponding mapping is no£ integrable. However a complexity reduction 
occurs and using the method developed in [1, 22] one finds for the generating function of the degree defined in Eq. 1 

f ^ = (1 - x){l - x - X 2 ) 

leading to a complexity 

A = 1 1.618034 
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4- Mapping of the six families of stable patterns of III A 4 



These results are verified by direct inspection. Using the methodology and notation of [22], the collineation C, the 
inverse /, the generic (i.e. arbitrary q) degree generating function / and therefore the complexity can be computed 
for Pi, Qi, P 2 and Q 2 ■ 

• for Pi and Qi the mapping are trivial. Pi leads to a linear mapping with generating function j^z, and Qi to a 
mapping (K^ = -1). 

1 1 \ f X \ ( X + (2 - q)xi 

1 1-qh M*i -{ 



Cp = 



C c 



1 1 

1 e„ 



-co 

Xl 



Xl 
-X 2 



• for Pi and Q 2 , one can also find explicitly the expression of the collineation and of the inverse. The various 
expressions depend of the parity of | (as mentioned q should be even) yielding four mappings all having the 
same degree generating function 

1 + x 



fix) 



(1 - x){\ - 2x) 



giving an integer complexity A = 2. Introducing e q = (—1) 2 the corresponding collineation and inverse are given 
below 

xl - (q - 2)x\ -{q- 4)x xi + e q x x 2 



Cp 2 = 




1 



-1 



/--even 
°Q 2 



- 1 ^ 
± 2 

1 

- 1 - 



2 1 
2 1 



c: 



odd 



1 1 



1 1 _ 

2 2 

1 

2 -1 + 1 _l 

2 2 ' ± ± 



1 



1 



reven 
J Q 2 



r odd 




The superscript refers to the parity of 



-(.t - e q x 2 )xi 
+ e q (q- 2)x\ + (q - A)x 2 xi - x x 2 



(2q - A)xl - 2qx\ + Ax x 2 
-4(x - x 2 )xi 
-(q - 2){q - 4)^ + q{q - 2)x\ - Ax 2 2 - i(q - 3)x X 2 

4(xi - x 2 )x 
2qxl + (4 - 2q)xj - Ax x x 2 
-q(q - 2)x 2 + (q- 2)(q - A)x\ + Ax\ + 4(g - 3)xix 2 
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IV. NUMERICAL RESULTS 
A. Computer-aided method 

To study the case where q is neither a prime nor the square of a prime, we use a computer. In principle there is no 
difficulty since it is "only" a matter of generating the patterns or the signed-patterns and check the stability. To test 
the stability we can use the formula Eq. 10-13. However, in practice, this would be tractable only for a very small 
value of q since the number of patterns grows extremely fast with q. 

If a partition £ = {Ei, ■ • • , E r } is product stable then for any k, Ek must be such that i, j € Ek => (Cy(x(-E , fc))) i = 
(Cy(x(Ek)))j- Since the cyclic matrices associated to a stable pattern form an algebra, we see that any subset Ek 
must be convenient (see Eq. 8). This simple remark tells that it is not necessary to generate all possible subset 
Ek, but only convenient sets: one can first enumerate the 2 q subsets E of {0, • • • , q — l}and keep only those which 
are convenient. The key point here is that a subset E is convenient irrespective of the way the remaining indices 
are grouped into others subsets. Then, using a tree structure, we can associate convenient subsets to make stable 
patterns. In order to minimize the tree structure, one can also use condition on pairs (or more) of index subsets, 
also deduced from Eq. 8. This procedure applies mutatis mutandis to inverse-stability when the matrices x{E) are 
invertible. Note that this procedure can also be applied to the search of non cyclic stable patterns. 

In the case of cyclic matrix we can go further and avoid the consideration of all possible subsets, retaining only 
the convenient sets. Indeed using the result of section II A 3 one can generate directly the admissible sets. We then 
consider these sets as "atoms" to be combined to produce the patterns. This can also be implemented using a tree 
structure. In the results shown below, we did not use this last remark. 
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q 


2 


3 


4 


5 


6 


7 


8 


# tested pattern 


3 


11 


49 


257 


1539 


10299 


75905 


P-Stable Pattern 


1 


2 


3 


3 


7 


4 


10 


P-Stable signed-Pattern 








6 





3 





17 


TP-Stable Pattern 








2 





3 





11 


IP-Stable signed-Pattern 








2 





1 





19 


Total 


1 


2 


13 


3 


14 


4 


57 



Table I: Number of P- and /P-stable patterns and signed-patterns. 



<7 


9 


10 


11 


12 


13 


14 


15 


16 


17 


# stable patterns 


7 


10 


4 


32 


6 


13 


21 


37 


5 


q 


18 


19 


20 


21 


22 


23 


24 


25 


26 


# stable patterns 


42 


6 


47 


28 


14 


5 


172 


13 


19 


q 


27 


28 


29 


30 


31 


37 


41 


43 


49 


# stable patterns 


25 


61 


7 


148 


8 


9 


8 


8 


21 



Table II: Number of product-stable patterns. 



B. Results 

In table I we present the number of P-stable and IP-stable patterns and signed-patterns for 2 < q < 8. In table 

II we present the number of P-stable patterns for larger values of q. All these numbers have been found using the 
algorithm presented above. We have verified that for q prime (i.e. q = 2, 3, 5, 7, 11, 13, 17, 19, 23, 29) it 
corresponds to the result of section IIIAl, and for q the square of a prime (q = 4, 9, 25) to the result of section 

III A 2. In Table I the first line is the number of convenient sets. The explicit expression of each pattern is not given 
in the text, but can be downloaded from the site [28]. However the case q = 8 is given in detail in Appendix A. 
Finally a maple program to generate all the stable patterns for q = p 2 with p prime can be download from [28]. 

V. CONCLUSION 

In this paper we have shown several results concerning stable patterns in the case of cyclic matrices. The notion of 
signed-pattern arises naturally when one studies IP-stability as a consequence of a duality between cyclic matrices 
and their Fourier transform. We find in particular an exact correspondence between IP— stable patterns and P— stable 
signed-patterns, which justifies, a posteriori, the introduction of signed-patterns in this cyclic matrix context. The 
main results Eq. 10-Eq. 13 enable to find all /-stable patterns and signed-patterns when the number of states is 
a prime or the square of a prime, and to find some, but not all, stable patterns for composite integer values of q. 
This provides examples of birational transformations of an arbitrary large number of variables. We have computed 
the complexity of the corresponding transformations in some cases, finding a complexity reduction. In particular we 
have recovered a one-parameter family of integrable transformations, for which we have given explicitly the matrix 
representation when it exists. The case of the monocolor /-stable signed-patterns has been solved, demonstrating a 
conjecture about Hadamard matrices in a particular case. We also present an algorithm to find /-stable patterns. 
Although this algorithm is exponential, it can be used for not too large values of q. 

It would be interesting to generalize our result to arbitrary value of q and to perform a more systematic analysis of 
the complexity of the associated birational transformation. Finally the same problem for non cyclic matrices should 
also be investigated, but it seems to us that it becomes a very complicated task, as a consequence of the loss of the 
discrete Fourier transform. 



Appendix A: STABLE PATTERNS FOR q = 8 

We list below the stable patterns for q = 8. The number before the eight letters between bracket is the arbitrary 
label of the pattern. The sequence of eight letters designates the first row of the cyclic matrix in the direct space, and 
the diagonal of the matrix in the Fourier space. When a letter is repeated (resp. negated) this means that the two 
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corresponding entries of the matrix are equal, (resp. opposite). For example the pattern number 10 below corresponds 
to 
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and the cyclic matrix associated to pattern 17 with first row [a, 6, a, b, c, 6, a, b] becomes by fourier transform the 

diagonal matrix M}° . to which one associates in the direct space the cyclic matrix Mj? , (see above). 
to tourier 1 J direct v ; 



Direct space — > Fourier space 
P-stable pattern — + P-stable pattern 

9[a, b, c, d, e, b, /, d] — » 9[a, b, c, d, e, 6, /, d] 
15[a, b, c, b, d, b, c, b] — > 15[a, b, c, b, d, b, c, b] 
22 [a, 6, 6, 6, 6, 6, 6, 6] -> 22[a, 6, b, b, 6, 6, 6, b] 
39 [a, 6, c, d, e, d, c, 6] — ■> 39[a, 6, c, d, e, d, c, 6] 
48 [a, b, c, d, e, f, g, h] -> 48 [a, 6, c, d, e, /, g, h] 
51 [a, 6, c, 6, d, e,c, e] — > 51[a, 6, c, 6, d, e,c, e] 
3[a, 6, c, 6, d, 6, e, 6] <-*■ 53[a, 6, c, d, e, 6, c, d] 
24 [a, 6, c, b, c, 6, c, 6] *-> 56[a, 6, 6, 6, c, 6, 6, 6] 



Direct space — > 
IP-stable pattern — 

4[a, a, a, a, 6, a, a, a] - 
6[a, 6, c, 6, a, d, e, d] - 
11 [a, b, a, c, a, d, a, e] 
14[a, 6, c, d, a, 6, e, d] 
17[a, 6, a, 6, c, b, a, 6] 
18[a, 6, c, 6, a, d, c, d] 
20[a, 6, c, d, a, d, c, 6] 
25[a, a, 6, a, a, a, c, a] 
29[a, 6, c, d, a, e, c, /] 
34[a, b, c, 6, a, b, d, b] - 
54[a, b, c, d, a, d, e, 6] 



Fourier space 
P-stable signed-pattern 

-> 55 [a, 6, —6, 6, —6, 6, —6, 6] 
+ 35[a, 6, c, d, e, — d, c, —6] 
— > 31 [a, b, c, d, e, —6, — c, — d] 
— ► 21[a, 6, c, —b, d, b, e, —6] 
-> 10[a, b, -b, b, c, b, -b, b] 
— ► 5 [a, 6, c, 6, d, —6, c, —6] 
— > 47[a, 6, c, —6, d, —6, c, 6] 
— ► 37[a, b, c, —6, — c, 6, c, —6] 
— > 13[a, 6, c, d, e, —6, /, — d] 
-> 49[a, 6, c, —6, d, 6, c, —6] 
— ► 44 [a, b, c, —6, d, e, c, — e] 



Direct space 
P-stable signed-pattern 

5[a, 6, c, 6, d, —6, c, —6] — 
10[a, 6, —6, 6, c, 6, —6, b] — 
13[a, 6, c, d, e, —6, /, — d] — > 
21 [a, 6, c, —6, d, 6, e, —6] — * 
31[a, 6, c, d, e, —6, — c, — d] — > 
35 [a, 6, c, d, e, — d, c, —6] - 
37[a, 6, c, —6, — c, b, c, — i>] — > 
44[a, b, c, —6, d, e, c, — e] — s 
47[a, b, c, —6, d, —6, c, 6] — 
49 [a, 6, c, —6, d, 6, c, —b] — 
55[a, 6, —6, b, —b, 6, —6, 6] - 



— ► Fourier space 
-> IP-stable pattern 
18[a, 6, c, 6, a, d, c, d] 
17[a, 6, a, 6, c, 6, a, 6] 
29[a, b, c, d, a, e, c, /] 
14[a, 6, c, d, a, 6, e, d] 
11 [a, b, a, c, a, d, a, e] 
■> 6[a, 6, c, 6, a, d, e, d] 
25 [a, a, 6, a, a, a, c, a] 
54[a, b, c, d, a, d, e, 6] 
20[a, 6, c, d, a, d, c, b] 
34[a, 6, c, 6, a, 6, d, 6] 
> 4[a, a, a, a, 6, a, a, a] 



Direct space — > Fourier space 

TP-stable signed-pattern — * TP-stable signed-pattern 

8[a, b, c, —6, a, —6, c, 6] — > 8[a, b, c, —b, a, —6, c, 6] 

16[a, —a, a, —a, 6, —a, a, —a] — > 16[a, —a, a, —a, 6, —a, a, —a] 

19[a, 6, c, d, —a, 6, — c, d] — ► 19[a, b, c, d, —a, 6, — c, d] 

23[a, 6, c, d, a, — d, e, —6] — ► 23[a, 6, c, d, a, — d, e, —6] 

27[a, b, c, 6, a, —b, c, —6] — * 27[a, b, c, 6, a, — b, c, —6] 

32[a, 6, c, d, —a, e, — c, /] — > 32[a, 6, c, d, —a, e, — c, /] 

38 [a, 6, c, d, a, —6, c, — d] — > 38[a, b, c, d, a, —b, c, — d] 

42[a, 5, c, —6, a, d, e, — d] — > 42[a, 6, c, —6, a, d, e, — d] 

50[a, 6, c, —6, a, 6, d, —6] — ► 50[a, 6, c, —6, a, 6, d, —6] 

l[a, b, c, 6, a, —6, — c, —6] <-> 52[a, 6, a, c, a, — c, a, —b] 

2[a, b, c, —6, —a, —6, c, 6] <-> 45[a, b, —a, c, a, c, —a, 6] 

7[a, 6, c, d, a, — d, c, —b] <-» 36[a, 6, c, b, a, —6, d, —6] 

12[a, b, c, —6, a, —6, — c, 6] <-> 41 [a, 6, a, —6, a, c, a, — c] 

26 [a, 6, —a, b, a, c, —a, c] <-» 57[a, 6, c, 6, —a, —6, c, —6] 

28 [a, 6, c, —6, a, —6, d, i>] <-> 33[a, 6, c, —6, a, d, c, — d] 

30[a, —a, 6, —a, a, —a, c, —a] «-> 40 [a, 6, —a, —6, c, 6, —a, —6] 

43[a, 6, —a, c, a, d, —a, e] <-> 46 [a, 6, c, d, —a, —6, e, — d] 
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In the following we show in detail that Eq. 13 is indeed verified on the example of the Fourier related pair of 
signed-patterns labeled 7 and 36 above. We used the letter E for pattern 7 and F for pattern 36, and u> 



2-TT 

exp =g-i. 







1 




V°7 



1 






\0J 



,x(£f) = 










w 





1 





1 

voy 






1 







,x(Ez) = 



fo\ 






1 



voy 



X(E+) 



f2\ 

2 

2 

2 

voy 



*x{Et 



( i \ 

u> 
i 

—Co 
-1 
— to 
— i 
V Q J 



,X{Ey 



( 1 \ 

Co 
—I 

—CO 
-1 
— LU 
I 

V w y 



/ 2 \ 


-2 

2 


-2 

V o y 



,x(£ 3 + ) 



/ i \ 

—Co 
— i 
u> 
-1 
Co 
i 

\-coJ 



/ 


1 






— w 












CJ 






-1 






w 






— i 




V 


—Co 


) 



X(E+ 



on another hand 



/2\ 


2 

2 

2 

voy 



X(F + ) 



, X (Et)-x(ET 



fl\ 




1 




voy 



/ o \ 

V2i 
2i 

V2i 


-V2i 
-2i 
\-V2i) 



,x{Et) = 



, X {F+)-x{F^) 



( \ 
1 


1 



-1 



v-iy 



/ 2 \ 


-2 

2 


-2 

V o y 



, X (E+)-x(Ez) 



,x{Ft) 



/o\ 



1 






voy 





V2i 
-2i 
V2i 


-V2i 
2i 

\-V2iJ 









1 

voy 



2 X (F+) + 2 X (F+) + 2 X (F+ 



It is then straightforward to verify that the four vectors of Eq. Al and Eq. A2 span the subspace: 
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Appendix B: STABLE PATTERNS FOR q = 25 = 5 2 

We list, as an illustration, the stable patterns for q = 25. There are 1 + r(4) + r 2 (4) = 1 + 3 + 3 2 = 13 such stable 
patterns. 

Firstly there is the simple pattern corresponding to the standard Potts model: 

1 [a, b, b, b, b, b, b, b, b, b, b, b, b, b, b, b, b, b, b, b, b, b, 6, b, b] 
then the pattern corresponding to the subgroup of Z25 = {1, 2, 3, 4, 6, 7, 8, 9, 11, 12, 13 • • • , 24} : 



L\ — Z24 

L 2 = {1,4,6,9,11,14,16,19,21,24} 

L 3 = {1,6,11,16,21} 

L 4 = {1,7,18,24} 

L s = {1,24} 

Le = {1} 



which gives the six patterns 



2 [a, b, b, b, b, c, b, b, b, b, c, b, b, b, b, c, b, b, b, b, c, b, 6, b, b] 

3 [a, 6, c, c, b, d, b, c, c, b, e, 6, c, c, 6, e, 6, c, c, 6, d, 6, c, c, 6] 

4 [a, 6, c, d, e, /, 6, c, d, e, g, 6, c, d, e, /i, 6, c, d, e, i, 6, c, d, e] 

5 [a, c, d, e, e, 6, /, c, /, 3, 6, d, 3, 3, d, 6, g, /, c, /, 6, e, e, d, c] 

6 [a, d, e, /, .9, b, h, i,j, k, c, I, m, m, I, c, k, j, i, h, b, g, /, e, d] 

7 [a, 5, /i, i, J, c, k, I, n, b, d, p, m, o, q, e, r, s, t, u, /, v, w, x, y] 



Finally the remaining patterns are computed fromZg = {1,2,3,4}: 



K 1 = {1,2,3,4} P~ 1 (K 1 )=Z* 25 

K 2 = {1, 4} P- l {K 2 ) = {1, 4, 6, 9, 11, 14, 16, 19, 21, 24} 
K 3 = {1} P~ 1 (K 3 ) = {1,6,11,16,21} 



yielding the six last patterns: 



8 [a, b, b, b, b, c, b, b, b, b, d, b, b, b, b, d, b, b, b, b, c, b, b, b, b] 

9 [a, 6, 6, 6, 6, c, 6, 6, 6, 6, d, 6, 6, 6, 6, e, 6, 6, 6, 6, /, b, b, b, b] 

10 [a, c, d, d, c, 6, c, d, d, c, 6, c, d, d, c, 6, c, d, d, c, 6, c, d, d, c] 

11 [a, 6, c, c, 6, d, 6, c, c, 6, e, 6, c, c, 6, /, 6, c, c, 6, 3, 6, c, c, 6] 

12 [a, c, d, e, /, 6, c, d, e, /, 6, c, d, e, /, 6, c, d, e, /, 6, c, d, e, /] 

13 [a, d, e, /, 5, 6, d, e, /, 5, c, d, e, /, 3, c, d, e, /, g, 6, d, e, /, 3] 
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